Condensed Matter Physics, ????, Vol. ?, No ?(??), pp.ffl-?? 



Grassmann techniques applied to classical spin 
systems 

M. CluseR J.-Y. FortirP, 

1 Department of Physics and Center for Soft Matter Research, New- York University, 4 
Washington place, New- York NY 10003, USA 

2 Institut Jean Lamour Departement de Physique de la Matiere et des Materiaux, Groupe de 
Physique Statistique CNRS - Nancy-Universite BP 70239 F-54506 Vandoeuvre les Nancy 
Cedex, France 

June 3, 2009 
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1. Introduction 

Classical an quantum spin models such as Ising model play an important role in 
the field of statistical physics as they allow for an accurate understanding of critical 
phenomena in general. Many techniques [ [U [2] were developed in order to deal 
with the difficulty of estimating the partition function and other thermodynamical 
properties in the critical region in dimension more than one. An exact mathematical 
description of the two-dimensional (2D) Ising model relies on the Jordan- Wigner 
transformation [[3] which maps the product of Boltzmann weights onto a fermionic 
action of free fermions with a mass vanishing at the second order critical temperature 
given in dimensionless units T c = 2/ ln(l+\/2) — 2. 2691851. Also a method based on 
the correspondence between the Ising model and dimer problems [H] uses the notion 
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of Pfaffians, which are directly connected to integrals over Grassmannian objects. 
Both fermions and Grassmann variables are therefore closely tied to the Ising model. 
A direct introduction of Grassmann variables as an alternative tool to solve the Ising 
model was done long ago in the 80's by Bugrij [|3] and Plechko [E] (see also a later 
discussion by Nojima [ 7j). It is based on a simple integral representation of the 
individual Boltzmann weights and which has the property to decouple the spins. 
The price to pay is a non-commutativity of terms arising from this representation. 
In order to deal with this particular representation, Bugrij used two families of 
Grassmann variables which commute with each other, then identified the resulting 
functional integral of the partition function with a determinant. From another point 
of view, Plechko introduced symmetries which order the non-commuting quantities 
so that the sum over the spins can be performed exactly. In this paper we review the 
process of how to generalize Plechko's method for Blume-Capel model [ 8], which 
is the simplest model beyond Ising, to spin-S Ising models and how to construct 
an exact fermionic action for each case. This would provide a natural extension of 
the exact fermionic quadratic action found for the spin S = 1/2 Ising model. In 
particular, we will build on previous work on the Blume-Capel (BC) case S = 1 [[8] 
where a line of second-order critical points is terminated by a tricritical point. This 
is the next case beyond the Ising model and which possesses a rich critical behavior. 
This model was used to qualitatively explain the phase transition in a mixture of 
He 3 -He 4 adsorbed on a 2D surface [[9]. Below a concentration of 67% in He 3 , the 
mixture undergoes a A transition and the two components separate through a first 
order phase transition with only He 4 superfluid. On a 2D lattice, He atoms are 
represented by a spin-like variable, according to the following rule: an He 3 atom is 
associated to the value 0, whereas a He 4 is represented by a classical Ising spin taking 
the values ±1. Within this framework, all the lattice sites are occupied either by an 
He 3 or He 4 atom. In addition to nearest-neighbor interactions, the energy includes 
a term A J2mn ^mni w ith = 0,1, to take into account a possible change in 
vacancies number. A can be viewed as a chemical potential for vacancies, or as a 
parameter of crystal field in a magnetic interpretation of atomic physics. It would 
be in particular interesting to have a fermionic description of the BC model in order 
to obtain more information about the kind of interaction fermions living on the 2D 
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lattice have compare to the Ising free fermion case. 

So one of the main question is how to obtain a generic fermionic action for a 
spin-S model and what does this method teaches us for the BC model in particular. 
We explain in the next section the main ideas of this method. 



2. Description of the fermionization for general classical spin-S 
models 

Let us consider the following Hamiltonian on a 2D lattice of size L x L 

L L L L 

H = — ^] J ^S mn S m+ i n + •S'mnS'mn+iJ + A 'S'mrn (1) 

m=l n=l m=l n=l 

where J is the Ising coupling constant and Ao the splitting crystal field or rep- 
resents a chemical potential in the BC case. In particular for Ao large and posi- 
tive, it favors small spin values. This crystal field can be replaced by any potential 
V(Smn) depending on the square of the local spin. Spins S mn take 2S + 1 val- 
ues with S mn = —S,—S + 1, • • • ,S. The partition function is the sum over all 
possible spin configurations Z = Tr exp(— (3H). Z contains products of the Boltz- 
mann weights exp(K S mn S m+ i n ) (where S mn and S rn+ i n are neighboring spins and 
K = J/k B T) which take q+1 = S(S + 1) + 1 different values if S is an integer, and 
q + 1 = (S + 1/2) (S + 3/2) values if S is half-integer. Since there are q + 1 possible 
values for each Boltzmann weight, we can project each of them onto a polynomial 
function of degree q in the variable S mn S m+ i n : 

q q 

+ %aSmnSm+ln) i (2) 

fc=0 a=l 

where the q+1 constants Uk are determined by solving the linear system of q + 1 
equations satisfied by the above relation. To see on specific examples how it works, 
let us consider first the Ising case, S — 1/2. Since S is half integer, we have q — 1. 
Therefore 



exp(KS mn S m+ln ) = ch(K/A) + 4 sh(K/A)S mn S m+ln , u Q = ch(X/4), Ul = 4sh(tf/4). (3) 
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In the Blume-Capel model, since S is integer, we have q = 2 and it is straightforward 
to show that 

exp(KS mn S m+ln ) = 1 + sh(K)S mn S m+ln + (ch(K) - l)(S mn S m+ln ) 2 , 

u = 1, m = sh(AT), u 2 = ch(K) - 1. (4) 

For S integer the first coefficient Uq is always unity, and from equation (j2J) we can 
write 

u k = u ^ x »i x a 2 ■ ■ -Xa k , ^-<k<q. (5) 

ai<ct2<--<a k 

We will set for convenience in the following Uk> q +i = since the polynoms are all 
finite. Our purpose is to transform the partition function Z which is a sum over spin 
variables into a multiple integral over Grassmann variables. For this let us introduce 
q pairs of Grassmann variables [El El [10] (a* n , a mn ) on each site for the horizontal 
direction and q additional pairs (& mn , & mn ) for the vertical direction. Here a takes the 
values 1, . . . q. There are therefore 4q Grassmann variables at each site of the lattice. 
In particular the Ising model is represented by two pairs of Grassmann variables per 
site which can afterward be reduced to only one pair [E]. For each couple of terms 

(1 -|- x a S mn S m ^ r \ n ^ (1 -|- x Q! 5' ?ri?1 >S , ?Tin -|_x) (6) 

appearing in the partition function, we introduce the following integral representa- 
tion 



1 _|_ T Q Q , . - / rln a dn a p a ™ a raJl _i_ n a Q Wl _i_ r 7, a <7 , , 

1 + x a S ' mn »S 'rnn+i = / db mn db mn e mn mn (l + 6^ in 5' mri )(l + x a b mn S mn ^i). (7) 



From the last expression, we introduce the link factors A^ n = 1 + a^ in S mn , A^ +ln = 
1 4~ x a CL mn S m +i n , B mn = 1 + b mn S mn , and B mn _^_^ = 1 + x a b mn S mn j r i, so that the 
partition function can be written as 



Z = i^ L2 Tr {5iajfe} 



J] e A5 ™ x [n(« +1 j n( 5 w 
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where A = — /3Aq. The mixed trace operator introduced in the last expression is 
defined by the following sums and integrals: 



Tr {5 ,a,fe}[.] = Tr {5} / [ JJ da^ n da^ n db^ n db^ n x e< 

mn,a 

Inside the integral symbols, the pairs of Grassmannian link factors in brackets 

(^mn^m+ln)) (^mn^mn+l) {") 

can be moved freely with the other terms since they correspond to commutative 
scalars after integration. In particular, it is convenient to rearrange the products 
over a in order to put together the link factors of different a with the same site 
indices (m, n) using the mirror ordering symmetry introduced in Plechko's method [ 
[6] in the context of the 2D Ising model, and which is still relevant in the spin-S case: 
i 

n( A a A a ) = (A 1 A 1 ) ( A q A q \ 
\ mn m+ln/ V mn m+ln) ' ' ' \ ran m+ln/' 



a=l 

{^■mni^-mn ' • ' (Aral (^mn^m+ln) ^m+ln) • •^m+ln) i 



a=l J \a=l 

where the arrows indicate that the product is ordered, i. e. increasing label a in the 
first product from left to the right and in the second one from right to the left. For 
convenience, we will use the notation 

Omn — W_A mw Om+in = J^J^4 m +i n , (9) 

a a 

for objects on the horizontal links and 

a a 

for the ones on vertical links. Then the partition function can be rewritten as 

z = u 2 l2 tt {s} j [n^^^^^^^^n^^ 5 ^ 1 ")^^ 

mn,a — 

= u 2 L2 Ti {s ^ b} Y[e AS ^(O mn O m+ln )(V mn P mn+1 



mn+1 , 
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At this stage of the algebra, we use the mirror and associative symmetries which 
were used for solving the Ising model [ El [TO] and which are still valid here to 
rearrange the operators O and V. In principle boundary terms should be treated 
separately in order to obtain the exact finite size partition function depending on 
boundary conditions [ 6] but they are not relevant in the thermodynamical limit 
L — > oo we are interested in here. Here we consider instead the simple case of free 
boundary conditions, and we obtain the exact expression after rearrangement of the 
O and V operators: 



Z = u 2 L "Tr { s,a,b} 



n(T\e AS ™(6 v o \T\v ) 
\ II ° I ^ran i mn\^mn III' ran J 

n=l m=l m=l 



Now, from this expression, the spins can individually be summed up from Sm to 
Si n for any given n. We will need to introduce the following weights W mn which 
include all the dependence on the individual spins S mn 



''mn / ° ^mn' mn^mni • n 

Smn — il 



= e e A5 ™n( i+c -M' (i2) 

S mn =±l a=l 

where we have defined the following Aq sets of Grassmann variables c n \ n in the 
following order: 

1 -q 2 _ -q— 1 q _ -1 

^mn ^ qQ'm—lni ^mn ^q—^m—lni ) ^mn -^l^m— In' 

r S+l _ r Ul r Q+2 = Zq-1 ... 2q _ 71 

c 2 9 +l = 1 2q+2 = 2 ... 3q = q 

3q+l _ Ll r 3q+2 = h 2 ... 4? _ w /i o\ 

mn mn' mn mn' ' mn ran' \ / 

The sum over SVnn = ±1 in equation f[T2j) can be performed by noticing that only 
products involving an even number of S" mn give a non-zero contribution. We also 
define the scalars (we remind that A = — /3Ao) 

s 

«*= E ^ n exp(A5L), (14) 
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and the ordered products 

Jk) _ „ai a 2 ,,. r a k JO) — i /i r\ 

Hmn / j ^mn^mn ^mni Hmn ) V / 

ai<a,2<—<a.k 

with qfnn = c x mn ■ • • c^ n the term of highest degree in Grassmann variables. Using 
these quantities, it is easy to show that the partial Boltzmann weights (|T2|) are given 
by 

2q 

W mn = J2 a KV™- (16) 

k=0 

Then the fermionic representation of the partition function can be expressed as a 
multiple integral over Grassmannian variables only 

Z = ul L2 Tv {Ml] \[W mn . (17) 

mn 

For small values of S, the weights W mn can be exponentiated so that a fermionic 
action can be defined. Indeed, since the first term of W mn is the pure scalar a and 
the others products of pure commutating Grassmannian objects, it is tempting to 
exponentiate the sum ffTol) to obtain directly a fermionic action. This comes from 
the simple observation that for any Grassmann variable a, we have 1 + a = e a . Of 
course, the exponentiation of the sum ffTBT) is more complicate. For example, for 
commuting objets a and b such as the q^s, we have 1 + a + b = exp(l + a + b — ab). 
In this case the order of the polynomial object inside the exponential is bigger than 
in the original sum since the extra counter-term ab is necessary for the identity 
to be exact. These weights are moreover connected by nearest- neighbor interactions 
hidden in the variables c^ n . In the case of the Ising model, where the exponentiation 
can be done quite easily, the argument of the exponential is purely quadratic in the 
c^ n 's and therefore the partition function can be integrated out with the use of 
a determinant or a Bogoliubov transformation in the Fourier space. Moreover, the 
Aq = 4 Grassmann variables in this case can be reduced to 2 by partial integration 
of non relevant variables. In the BC model, the argument is a polynomial of degree 8 
in Grassmann variables since there are 8 independent variables (4g = 8). In general 
we expect naturally the argument to be at most a polynomial of degree 4q in these 



7 



MC-JYF 



variables, which can be reduced or not by partial integrations. Except for the case 
q — 1 however the partition function can not be expressed as a determinant, so that 
a full exact solution of the partition function can not be found this way. If the action 
is quadratic, the use of the following Gaussian integral [ d2] , defined on Grassmann 
set of variables {a i? <ii}i=i,..,N, and for a square matrix A 



allows us to express the partition function as a determinant. Quadratic fermionic 
form in the exponential (fT8l) is typically called action for a free-field theory. When 
the action is non-quadratic, the integral is not Gaussian and can not be expressed 
as a determinant, which yields in principle to a non integrable theory. However, 
physical information such as bare masses (see last section) can be extracted from 
these non-quadratic actions which represent generic theories of interacting fermions. 



3. Fermionic action of the Blume-Capel model 

In this section, we consider the case S = 1 (Blume-Capel model) which is the 
simplest example of a classical spin beyond the Ising model. It possesses in the 
phase diagram (T, Ao) a second-order critical line separating a ordered phase from 
a disordered one and terminated by a tricritical point (see figured] below). From the 
previous section equation (jTTJ) allows us to write an action after exponentiation of 
the Grassmann variables which can be done exactly after some tedious algebra. The 
4 pairs of variables per site can however be reduced to 2 pairs by partial integration. 
Another simpler way of obtaining this BC fermionic action is possible [[11] using the 
Z2 symmetry of the spin variables S mn . Indeed the partition function is invariant if 
we perform the gauge transformation S mn — > o mn S mn with a mn = ±1. In this case 
it is possible to simplify the process of the previous method and write an action 
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containing only 2 pairs of variables per site instead of 4: 

/L L L L 

II II da mn da mn db mn db mn exp j ^ ^ (19) 
m=l n=l m=l n=l 

^•mn^mn b m nbmn Q"mnbmn ^(^m— In bmn—ljiflmn b m n) t dm—lnbmn—l 
90 Q'mnO'mnbmnbmn exp ( '")'Cl J71 _i n G(, Jn _i n 'y&rren — l^mrj. — 1 t Q"m— Inbmn— l) 

where we have introduced the following constants: 



1 

2 cosh 2 if' ' ~ cosh if 



9o = 7 72 — > 7=1- — - = 1 - Vl ~t 2 , t = tanhif. (20) 



The fermionic integral (Tl9j) is the exact expression even for a finite lattice, provided 
we assume free boundary conditions for both spins and fermions. The other possible 
form for the partition function with periodic boundary conditions in both direction 
can be written in a similar way as the Ising model on a torus [El [T3l [Hj. The par- 
tition function would be the sum of 4 fermionic integrals with periodic-antiperiodic 
boundary conditions for the fermions. In the expression ( fT9l) . we can recognize the 
sum of the Ising action, which here appears as the Gaussian part of the total action 
[El EH: 

L 

•Slsing ^ ^ &mn®mn b mn b mn -\- d mn b mn -\- t(fl m _i n -|- &mn-l) ("-mn b m n) t ^m— ln^mn— 1 ; 
m,n=l 

and a non-quadratic interaction part, which is a polynomial of degree 8 in Grassmann 
variables (which can be seen if we expand the exponential inside the action): 

L 

S\nt 90 ^ ^ "mn'^mn^mii^mn ^Xp ^ T^m— ln^m— In 7^mn~l^mn-l t O-m— ln^mn — 1^ -(21) 
m,n=l 

This allows us to rewrite the partition function as a fermionic field-theory in a 
compact form 



Z = (2e A cosh 2 K) l2 J VaDaDbVb exp(Si sing + S int ) . 



(22) 



The BC model differs from the Ising model by the interaction term in the action 
(1211) which is not quadratic. Therefore the BC model is not solvable in the sense of 
free fermions as a determinant of some matrix, unlike the 2D Ising model. 
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3.1. Mixed representation of the BC model 

The coupling of Grassmann variables in equation (|2"T|) prevents us to integrate 
further and reduce the number of variables per site unlike the Ising model where the 
minimal action contains one pair only [[18j [7] . The minimal action of the Ising model 
admits an interpretation in term of Dirac representation of free fermions which be- 
come massless at the critical point. In a previous work we were able to reduce the 
number of Grassmann variables by partially introducing hard core bosons in the 
previous action, since terms such as r] mn = a mn a mn or r mn = b mn b mn may have an 
interpretation of local densities or occupation numbers. Variables 7] mn and r mn are 
commuting and nilpotent, 7^ n = r mn = 0. We can replace the quantities depending 
on a mn a mn and b mn b mn , especially in the interaction part, by their respective nilpo- 
tent variables, using, for this task, a general definition of Dirac distribution for any 
polynomial function / of a mn a mn or b mn b mn [HI]: 



A natural definition [US] of the integrals involving commuting nilpotent variables is 
to impose the following rules (and similar for fj mn , f mn ): 



This change of variables allows us now to integrate over the a mn 's and 6 mn 's in the 
new action. One advantage is that after this operation there are only two fermionic 
variables per site, although two additional pairs of bosonic variables have been intro- 
duced. In fact we can integrate over one pair of bosonic variables [ITTj. for example 
Vmn, fmn, using the help of integration rules and Dirac function given by (1231) . At 
the end, it remains a mixed action made of one pair per site of fermionic and bosonic 
variables respectively, with an interaction between fermions and bosons. A conve- 
nient replacement of the variables a mn by c mn and b mn by — c mn in the final integral 
leads us to isolate the minimal local action for the pure Ising model [ [161 E] with 
one pair of Grassmann variables per site: 




(23) 




(24) 




(25) 
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and the interaction part 



Sint — go ^mn^mn [(1 - T 7 7m-ln)(l - 7 r mn-l) + t 2 C m -i n C mn -i] , (26) 

with the quantities 

Qmn C rnn C mn -|- tc mn {c m —\ n 

Qmn C mn C mn ~\- tC mn (^C m ~\ n 

The Ising part is the same action that results from the integration over a mn ,b mn 
from the original Ising case. The introduction of nilpotent variables was necessary to 
achieve this partial extraction of the Ising contribution. The physical interpretation 
of the previous mixed representation is that it can be possible to describe the BC 
model with fermionic variables for the states S = ±1 and bosonic ones for states 
S — 0. In the limit A — > — oo, the system is completely described in terms of 
fermions (Ising sector), while when A is increasing fermions and bosons begin to 
interact. Beyond a critical value of A , fermions form bosonic pairs and in the limit 
A — > +00, all fermions condense into bosons, leading to a purely bosonic system. 
This view should be supported by further analysis. 



Cmn-l) C mn \c mn -|- t(c ?n _i n C mn _x)] , 
Cmn-l) [Cjnri ^(Cm— In Orm— l)]Cmn . (27) 



3.2. Corrections to the effective action in the continuum limit 

The integration of the previous action fT26l) over variables (r] mn , r mn ) can be per- 
formed perturbatively, as part of an expansion in the low momentum limit. We will 
define formally the derivatives of Grassmann variables [ITS], d x c mn = c mn — c m _i„ 
and d y c mn = c mn — c mn _i in the limit of large L. In this limit and in the Fourier 
space, the high order derivatives account in the action for a small contribution in 
momenta k = 2ir(m, n)/L, with m, n -C L positive integers. We would like to obtain 
in this limit the non trivial part of the non-quadratic interaction in term of variables 
Cmn, Cmn only. The procedure is described in reference [ [TTJ and based partially on 
substitution rules such as 



TjmnTmn * C mn C mn , f] mn > Q m n j T~mn * Qmn • (28) 
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There are unfortunately more complicate terms in the resulting action than by using 
the substitution rules alone, such as 

but they can be discarded in the approximation scheme above in the sense they 
correspond to corrective terms higher than quartic polynomials or quantities of the 
order of 0{qq) where go is the natural parameter of the expansion. It is exponentially 
small in the region where Ao is large and negative (Ising behavior). At the lowest 
order we found that the effective action (I25II26P can be approximated by the following 
expansion with respect with g 

^effective Rising -|- (jo ^ ^ C mn C mn |^(1 7?m— ln){l ""YQmn—l) t C m _i n C mJ1 _iJ 
m,n 

~^~9o^y ^ > -\- . . . . (30) 

m,n 

From the previous result, it appears to be suitable to express the quadratic and 
quartic parts in the Fourier space (in the large but finite L limit), where we define 
the following transformations 

c(r) = — ^ c k exp(ik.r) , c(r) = ^ c k exp(-zk.r) . (31) 

k k 

The Ising part of the action can be written as 

Rising = ^t mBC + ^ + l ^ kx ~ ^)]( c kCk - c_ k c_ k ) + 2itk x c ]i c-y i + 2#fcj,C k C_ k ,(32) 

keS 

with wbc = 1 — It — t 2 + go and the quartic term can be express as 



Sint = 90^2 ^ / k",k"-qCk'Ck"Ck'+qCk"-q, (33) 

with the potential 



L 2 

k',k",q 



Vk,k' = -ak x k' y + a'(k x k' x + k y k' y ), 

a = t(t + 2 7 ) , a' = 7(1 - t) . (34) 
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We notice that the bare mass of the theory is given by 



m B c = 1 - 2t - t 2 + g = m Ising + g , (35) 

where mi S i ng = 1 — It — t 2 is the Ising mass which vanishes at the critical value 
tanh(fT c ) = \/2— 1 corresponding to the second order transition point T c = 2/ ln(l + 
\/2) ~ 2.2691851 in units of J/ks- In the BC model, the critical temperature is 
shifted by the parameter g Q which depends on the temperature and A . The lo- 
cation of the second order critical line goes from the previous Ising critical value 
T c = 2/ log(l + v^2) when A — > — oo to the zero temperature point (T c = 0, A = 2) 
continuously where the transition can be proved to be first order by a simple ener- 
getic argument. In figure ([T]), we have reported the critical line given by mec = 
and the different numerical results found in literature [ITT]. In general the agree- 
ment is good, which validates the fermionic theory giving a bare mass vanishing at 
locations close to critical point values found by numerical methods. The presence of 
a tricritical point is induced by the interaction term (1341) which renders the second 
order line instable. To see why, let us consider the infrared limit on the critical line. 
The spectrum is given by the lowest terms of an expansion of the effective action 
with respect with kinetic terms, an the contribution to the partition function, in 
the Fourier space, is the product of partial integrals Z k such as Z = rj k Z k , up to 
the second order in the momentum k. For the Ising model and for small momenta, 
the factors are exactly of the form (mi sing + Ak 2 ), with A a constant equal to 
t(l — t 2 ) [HE]- The coefficient in front of the term k 2 in factors Z^ can by described 
physically as a stiffness coefficient. For the Ising model, the stiffness is always strictly 
positive even at the critical point. In the BC case, however, we have a line of critical 
points as Ao varies from negative to positive values up to Ao = 2. The effect of the 
interaction potential fl34|) is to modify the expression of the stiffness, which now is 
no more constant but depends on the angle of the vector k and also the temperature 
and Ao (see reference [ [TTJ for explicit details). The result is that in the BC case the 
effective stiffness coefficient vanishes at some point on the critical line, at a value 
close to Ao = 2, which indicates eventually the presence of a tricritical point. It can 
be shown that the partition function can be indeed written as a product over the 
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Fourier modes Z = likes ^ k w hh 



Z k = ml c + k 2 [A + B sin 29 k ) , 



(36) 



6k being the angle of the vector k, and A and B depending on temperature and Ao- 
As long as \A\ is larger than \B\ on the critical line, the transition is second order. A 
singular point can be reached if A 2 = B 2 , in such case Z k are not all strictly positive 
if mec = 0. Beyond this point the effective action (130]) is not sufficient to describe the 
critical properties of the model. If we compare the fermionic description of the BC 
model to a bosonic Ginsburg-Landau $ 6 theory describing first order transitions, the 
presence of a tricritical point would be equivalent to the fact that both coefficients 
of $ 2 and $ 4 terms vanish. 

3.3. Critical behavior of the BC model: diagrammatic expansion 

In this section, we further analyze the influence of the interaction potential V^k' 
on the renormalized mass, in particular the shift of the critical temperature which 
was in reference [ [TT] assumed to be given by the point where the bare mass mec 
vanishes. We would like in particular to apply diagrammatic expansion of the effec- 
tive action (I3T)1) . For this, it is useful to express the Ising part of the action in term 
of Nambu-Gorkov representation of the fermions [ [221 121] , using the two-component 
objects 



Formally, the Green functions can be defined within this representation by 2 x 2 
matrices 




(37) 



G(k) = ($ k $ k r 3 ) 



(38) 



where t% is the Pauli matrix 




(39) 
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Figure 1. (color online) Figure taken from reference [ [TT] showing comparison 
between critical line denned by the vanishing mass (|35p (plain red line) and nu- 
merical results from Monte Carlo simulations. The black filled dots are from figure 
1, da Silva et al. [[19] (Wang-Landau method). The cross symbol indicates the 
tricritical point identified by [ [TT] using a Hartree-Fock-Bogoliubov approxima- 
tion for the quartic part (|33p of the effective action. The blue diamond symbols 
are from Ref. [[20]) the magenta triangles from Ref. [[H], and the green squares 
from Ref. [[22] (see also Table 1 for other numerical values at Aq = 0). 



The unperturbed part of the Green function Go is evaluted using the elements of 
the non diagonal but quadratic Ising action 



— 1 / nik 2itk y 

Go(k) = W + «»*,*, \2itk x m t 1 ' <40) 

where the momentum-dependent mass is defined by = itl-qq + it(t + l)(k x — k y ). 
The inverse is given by 



G -' ( k) = | ~ m ; 2ltk " i . (4i) 
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k" - q k' + q 




Figure 2. Representation of the interaction part of the potential. The blob rep- 
resents the potential interaction with incoming vector k" and outgoing k" — q. 

With this representation and the unperturbed Green function, we can write the Ising 
part as 

Rising = J2 ^T 3 G \k)<S> k , (42) 

kes 

where the set S contains half the momenta of the Brillouin zone. It is defined by the 
rule that if k G S, then — k does not belong to S. The interaction part can be put, 
after some algebra, into the following form 

Sint = -90^ Yl (^k"-qT 3 Vk"-q,k"^k") (*k'+qT3$k'), (43) 
k',k",q 

where the sum is not restricted to the ensemble S. We define the potential matrix 
V by 

Vic"-q,k" = • 44 

\ y k "_ q , k - J 

The two diagonal elements of this matrix are not equal since Vk,k' is not symmetric 
by exchange of the two momenta k and k' except when k = k'. 

We now perform a diagrammatic expansion with respect with g of the perturbed 
Green function G(k) which will allow us to compute the corrections to the mass, 
i.e. the shift of the critical temperature, by identification of the diagonal elements of 
the inverse propagator f(k) = G(k) _1 . The graphical representation of the matrix 
potential in term of diagram is displayed in figure [2j To do so we follow formally the 
Feynman rules which lead to the Dyson equation of the inverse-propagator in term 
of the self-energy S(k): 
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k k 




Figure 3. The four diagrams appearing at the lowest order in go- Only diagram 
(b) contributes to the mass in the low momentum limit k — ► 0. 



r(k) = G(k)- 1 = Go(k)- 1 - E(k). 



(45) 



The first terms contributing to the self-energy are given in figure [3J 



S(k) = -^^^[WqGolk + ^ + Golk-q)^-^ 



-0o 



7^ [^,kTrG (q) + f q>q TrG (q) 



(46) 



The renormalized mass m# is given in the limit when k is zero by the diago- 
nal components of the inverse-propagator rn(0) = r 2 2(0). In this limit, only one 
diagram is not vanishing, which corresponds to the diagram (b) of figure ([3]): 



i 

m R = T u (0) =m BC + go—^V^qTrGoiq), 



m BC - m BC x g 



AL 2 

—y 

2L 2 ^ m 2 BC + t 2 (l + t) 2 (q x - q y ) 2 + At 2 q x q y 



q,q 



• (47) 



The last sum over q can be evaluated in the continuous limit L — > oo. Setting 
q = 27r(^, we define for large L the two following integrals 

-27T /-27T 



h{m BC ) 
h(m BC ) 



go 



in' 



go 



JO 

2?r /-27T 



dq x dq y ^ 



QxQy 



m 



2 BC + t 2 (l + t) 2 (q x -q y ) 2 + At 2 q x q y 



o jo 



dq x dq 



ll + Qy 



y m 2 BC + t 2 (l + t) 2 (q x - q y ) 2 + At 2 q x q y 



(4f 
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These two quantities are finite when mec vanishes. To see why, we can consider 
polar coordinates q x = qcosO and q y = qsinO, so that, near the origin q = the 
second integral for example behaves like 

M-bc) ex ^ H ^ {mBc/(lt{1 + W + 1 + + 2/(1 + tf) M20) ^ 

When m-QQ = 0, this integral is finite since there is no singularity in the denominator. 
Indeed when 9 = ±7r/4 the modulus of the term — 1 + 2/(1 + 1) 2 is strictly less than 
one on the critical line. This would not be the case if the last term 4t 2 q x q y in 
the denominator and coming from the off-diagonal part of the Green function was 
absent. In this case the integrals would be singular in the limit of small mass mBc, 
the denominator would instead be equal to {mBc/qt(l + t)} 2 + 1 — sin(2#), and 
the singular part would behave like l/|m,Bc| by a simple scaling argument, which 
would cancel the other mass term mec in front of the integrals (H7|) . Then the 
renormalized mass would be shifted by a finite quantity, as well as the critical line. 
Here the renormalization only concerns the total coefficient of mec and this does 
not affect the critical line location: 

m B ~m B c(l-a/i(0) + c/I 2 (0)). (50) 

A plot of the positive ratio m^/mec evaluated at criticality mec — as function 
of A is given in figure HI It is close to unity for almost all values of A . 

This analysis shows that corrections to critical temperatures are indeed small 
and are of order higher than g . A further analysis should be carry out to obtain a 
finite shift by considering more complex diagrams. It is also supported by analytical 
values of m B c = which are close to numerical results reported in figure [TJ 

4. Extension to other spin-S models: generalization of the bare 
mass 

The previous bare mass computed in the BC model, mec, allowed us to obtain 
a precise description of the second-order critical line in the phase diagram. This was 
obtained by taking the limit of low momentum in the effective action, see (1301) . This 
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Figure 4. Plot of the coefficient renormalizing the BC bare mass as function of 
Ao on the critical line. The values are close to unity except near the region of the 
terminating point Ao = 2 where qq is not small. 

can be generalized for any value of the spin S, in particular for higher values of S. 
The equations obtained in section [2J equations (TToT) and (TT7I) . are general for they 
represent the fermionization of any spins-S model. 

The construction of the fermionic action is however not an easy task, unlike the 
BC model which is a simpler case, but we expect to be able to extract a bare mass 
associated to non kinetic terms, or term involving derivatives with respect with space 
variables. At first approximation, we assume that the partition function and the free 
energy are singular in the low momentum limit when this bare mass vanishes. In the 
continuum limit, the c's coefficients defined by relations (fT3|) can be rewritten using 
formal derivatives, such as c l mn = x g (a^ n — c^a^J, etc... The derivatives contribute 
only to the kinetic energy and not to the bare mass. Keeping the first terms of the 
expansion, c^ nn ~ x q ofl mn (as well as for the other coefficients c's), the weights W mn 
become uncoupled in the sense they contain variables depending only on local site 
(m, n) and we define the mass mg as: 



with ci,„ ~ x q a q ■ ■ ■ ^ 



CL„ — x\a mn and 



r h q 



r 1q 



x x b l mn . The 
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integral can be evaluated exactly by noticing for example that the arguments of the 
exponential b^ n J)^ nn can be combined with a a^ in (a; a a^ in ) that appears in some of the 
Qic=l 2q products to give a contribution x a . Indeed using the Grassmann integration 
rules f da. a = 1 and f da.l = 0, we can write 

/dn a dn a dh a dh a pCC+CCn" ( T 7, a \ - 

/ dd a da a db a db a ( 1 + a a a a 4- b a b a 4- a a a a b a b a \a a (r a a ) = 

/ uu mn uu, mn u,u mn u ' u m,n\ 1 ~ u mn Uj mn ~ u mn u mn ~ u mn u, mn u mn u mn I u, mn\- h aU' mn ) 

/ dd a da a db a db a b a b a a a (r a a ) = r 

Also a term a^ n a^ n can be combined with 6^ in (x Q &^ in ) to give the same contribution. 
Since the are ordered, there are also signs to take into account and coming from 
the fact the variables c^ nn have to be moved in the correct order before integration. 
We obtain after some algebra the general relation 

2q 

m s = ^2a k R k , (51) 

k=0 

where we have define the following quantities with initial condition R = Uq, 

k 

R k = y~"; uiu k -ia(l, k-l), (52) 

1=0 

and a(k,l) — 1 if k and I are both even, and a(k,l) = —1 otherwise. We can 
apply this result to different cases to check the validity of this relation. For the 
Ising model (S = 1/2, q — 1) u — ch(K/4) and u\ = 4sh(fT/4), we obtain mi/ 2 = 
2cosh(A/4)(w2 -uqUx/2 -uf/16), or 

m 1/2 = 2e A / 4 [l-sh(iT/2)], (53) 

which vanishes at the Ising critical temperature T c ~ 0.567296 or with the normal- 
ization t c = T c /S 2 = 2.269185, which is independent, as expected, of A . For the 
Blume-Capel model (S = 1, q = 2) we have mi = 1 + 2 exp(A)(l — 2ui + 2u 2 — u\ — 
2u\Ui + tt|), or more explicitly 

mi = 1 + 2e A [l - sh(2AT)]. (54) 
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Table 1. Comparison of critical temperature values at Ao = 0, solutions of equa- 
tion (I57p . with other methods found in bibliography. 



Spin S 


O I/O 

o = 1/2 


C 1 

o = 1 




O = o/z 


o = 2 


9 


1 


2 




rr 



r- 




T c 


0.567 296 


1.673 971 




3.277 561 


5.351 248 


4- T 1 1 C2 
tc = Ic/O 


2.269 loo 


1.673 971 




1.456 694 


1 O OT O 1 O 

1.337 812 


rveis. 


9 9RQ r PT^l 

z.zoy [ zoj 


i.ooy [ izoj . i.oyo 

1.694 [ES], L681 


r n-Ql 
[E7] 


1.401 [ IZ 1 1 \Zol IZoj 


l.OOO [ IZO] 


Spin S 


5 = 5/2 


5 = 3 




5 = oo 




q 


11 


12 




oo 




T c 


7.890 888 


10.894 806 




oo 




t c = T c /S 2 


1.262 542 


1.210 534 




0.925 148 




Refs. 


1.257 [[25] 


1.203 [[25] 




0.915 [[251 [29] 





This mass is directly proportionnal to the mass m-sc found in the previous section. 
Indeed, we have the relation 

mBC = gomi (55) 

and therefore both masses vanish on the same line of critical points. The coefficient 
<7o comes from a global rescaling of the Grassmann variables in the original weights 
W mn which leads to the coefficient g^ 1,2 in the BC function partition (jl~9"I) and f[22|) . 
instead of the coefficient {uq) 2l2 = 1 in front of f[T7|) . For Ao = we find in particular 
that t c = 2/arcsinh(3/2) ~ 1.673 971 (see table [J). 

The other masses are deduced by iteration of formula (loll and (|52|) . For higher 
values of spin S, we find: 

m 3/2 = 2e A / 4 [l - ak(K/2)] + 2e 9A/4 [l - sh(9#/2)] + 2e A/4 [l - sh(tf/2)], 

m 2 = l + 2e A [l-sh(2 J ftr)] + 2e 4A [l-sh(8 J fs:)], 

m 5/2 = 2e 25A/4 [l - sh(25AT/2)] + 2e 9A/4 [l - sh(9AT/2)] + 2e A / 4 [l - sh(K/2)], 

m 3 = l + 2e 9A [l-sh(18ir)] + 2e 4A [l-sh(8ir)] + 2e A [l-sh(2ir)]. (56) 
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For general spin S, we can extend the previous results to the formula 

s 

m s = e Au2 [l-sh(2a 2 K)}. (57) 

a=-S 

This is a simple result giving a precise location of the second-order critical lines by 
solving the equation of the bare mass = 0. We give for comparison tabulated 
values of t c at A = in tabled] for different values of S, and references to numerical 
results (Monte-Carlo simulations, high-temperature expansions) given in the liter- 
ature. In general, the agreement is good up to 1% in most cases. For half-integer 
values of S, the models possess an asymptote in the (T/S 2 , A ) plane. Indeed, for 
S = 3/2 for example, the equation given in f[5T)j) predicts the solution 

9t, r l-sh(2/9*)i 
A = - ¥ log[- r -^ ], (58) 

which is bounded by t c = 2/91og(l + y/2) = 0.252131 below which there is no 
continuation of the second-order critical line. In the limit of large integer S, the model 
defined in equation ([!]) is described by a continuous variable — 1 < x mn = S mn / S < 1 
and is called continuous Ising model. We can still obtain a finite value of the critical 
line by taking the asymptotic value of equation ( 1571) 



fv 1 r i 

m s>>1 ~ SV2t / dx e ~ Aox2/2 1 - sh(x 2 ) , (59) 
Jo L - 1 

and in particular for A = 0, we have the following expansion for large S 

m s (t, A = 0) ~ a(t)S - 1 + A + _1_ + . . . (60) 

with a(t) = 2 — 2-^/2/t sh(x 2 ) dx. We observe that the rescaled mass ms/S 

vanishes in this case when t c = 0.925 148, in good agreement with numerical works 
for this model [ [251 129] , and it is worth noting that equation (|59|) also possesses a 
non trivial solution at t — which is simply given by A = 4/V3 = 2.309 401. This 
value is different from the value 2 expected for all finite S models [ [8] . It can be 
suggested that there also exists a tricritical point before this non physical value is 
reached. 
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5. Conclusion 



In this review paper, we have presented a method which tries to operate a corre- 
spondence between classical spin models and fermionic systems. We have extended 
Plecho's method [ El [HI] based on Ising model to generalized spin-S systems. The 
method is based on the projection onto q polynomial components, q depending 
specifically on the value of spin S, of Boltzmann local weights given by equation (j2J). 
Then the introduction of 2q pairs of Grassmann variables per site and the use of spe- 
cial symmetries such as mirror and associative symmetries in 2D for Grassmannian 
objects allows us to perform exactly the sum over the spin variables. This gives a 
representation of spin-S models in term of fermionic multiple integrals (TTTj) . Effective 
actions can in principle be deduced from this representation. We have shown that 
such action can be built exactly for the Blume Capel model S = 1 (1301) and the bare 
mass (j55|) gives accurate description of the second-order critical line. We have seen 
that there is no shift of this mass due to the effect of quartic potential of the effective 
theory at the lowest order expansion in the coupling parameter go, implying that 
corrections to the critical temperature may be indeed small. This quartic potential 
is however responsible for the presence of a trier itical point, rendering the second 
order line instable by changing the sign of the stiffness coefficient or making the 
free-fermion spectrum itself instable. For general spin-S model, the bare mass can 
also be generalized and calculated directly in the low momentum limit fl51~l) without 
knowing the full effective action, and still gives accurate description of second-order 
critical points even in the limit of the continuous Ising model. 
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